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Reactive cross-sections for the collision of open shell S('D) atoms with ortho- and para- hydrogen, in the ki- 
netic energy range 1-120 K, have been calculated using the hyperspherical quantum reactive scattering method 
developed by Launay et al. [Chem. Phys. Lett. 169, 473 (1990)]. Short-range interactions, described using 
the ab initio potential energy surface by Ho et al., were complemented with an accurate description of the long- 
range interactions, where the main electrostatic (~ R~ 5 ) and dispersion (~ R~ 6 ) contributions were considered. 
This allows the comparison with recent experimental measurements of rate constants and excitation functions 
for the title reaction at low temperatures [Berteloite et al, accepted in Phys. Rev. Lett., 2010]. The agreement 
is fairly good. The behavior in the considered energy range can be understood on the average in terms of a 
classical Langevin (capture) model, where the centrifugal barriers determine the amount of reactive flux which 
reaches the barrierless transition state. Additionally, the structure of the van der Waals well provides temporal 
trapping at short distances thus allowing the system to find its way to the reaction at some classically-forbidden 
energies. Interestingly, the cross-section for para-hydrogen shows clearly oscillating features associated to the 
opening of new partial waves and to shape resonances which may be amenable to experimental detection. 



I. INTRODUCTION 

Until very recently, most experimental and theoretical stud- 
ies in reaction dynamics had focused on thermal (or higher) 
collision energies. With the field of cold molecules emerg- 
ing as an active interdisciplinary area of research, a great 
interest in understanding cold and ultracold chemical reac- 
tions has risen [lLl2l]- Electric and magnetic decelerators, he- 
lium buffer gas cooling, and magneto- and photo-association 
(among other techniques) are opening the access to cold and 
ultracold samples of molecules JH]. Given the success in 
controlling interactions achieved during the last two decades 
in the field of ultracold atomic physics j4], the question is 
whether a similar level of control is possible for more gen- 
eral processes involving molecules. Very few partial waves 
contribute at low temperatures. However, reactive processes 
have been shown to occur rapidly being dominated by tunnel- 
ing and Feshbach resonances, in a regime where the smallest 
barrier to reaction surely exceeds the collision energy JH-0]. 
Remarkably, at low kinetic energies the presence of external 
fields modifies the dynamics in an essential way, appearing 
as "knobs" to steer the system to the desired outcome llfJl . In 
addition, the theoretical advances in the coherent control of bi- 
molecular processes achieved in the nineties lfTll IT2I1 found in 
the thermal averaging of the colliding partners a fundamental 
obstacle. Consequently, the forthcoming access to coherent 
samples of molecules at very low temperatures has renewed 
the quest for the control of chemical reactions |[l3ll . 

Given the success in the production of (ultra-)cold alkali 
metal dimers, attention has mainly focused in such particular 
systems Ill4l4l9ll . In a landmark experiment, reaction rates for 
collisions involving fermionic 40 K 87 Rb, at temperatures less 
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than 1 /jK, were recently measured [20]. In contrast, advances 
in the analysis of other type of reactions are still hampered by 
the lack of versatile methods to produce molecules at low tem- 
peratures and the very low densities achieved. Luckily, emerg- 
ing experimental techniques like the Crossed Stark decelera- 
tors 12111 promise to provide detailed information on cold bi- 
molecular collisions. More related with our work, CRESU 
(Reaction Kinetics in Uniform Supersonic Flow) experiments 
implemented in Rennes B221I23I1 . and crossed molecular beam 
experiments performed in Bordeaux lf24l l25ll . are exploring 
barrierless reactions between neutral species in conditions 
aproaching the cold regime ll26ll . In particular, rate constants 
and cross-sections have been recently obtained for the title re- 
action at temperatures and kinetic energies as low as ~ 5K 
(1K~ 8.3 x 10" 3 KJ moL 1 ) JUl. We will analyze their results 
below. 

Given the availability of experimental results, the title reac- 
tion provides us with a good chance to check and extend our 
predictive power from the explored thermal regime, where ab 
initio reaction dynamics is routinely applied, to the low en- 
ergy one where such methodology should fail due to the lack 
of accuracy. In the limit, the deep ultracold regime, governed 
by Wigner laws, has been found amenable to parametrization 
in terms of a very few inputs (the scattering length, for ex- 
ample) which the experiments may provide soon 1 20] . Para- 
doxically, the extreme sensitivity of such parameters to small 
details of the PES, and the action of the surface as a whole, 
makes extremely difficult to predict them or, conversely, to 
deconvolute from their measurement the underlying interac- 
tions. Collisions in the range of ~ 1 K thus lie in the limits 
of what can be done using our conventional theoretical tools. 
Besides, such collisions are more than an obliged milestone 
in our route from the thermal to the ultracold regime, hav- 
ing interest on their own. They report the combined game 
of short-range (SR) and long-range (LR) interactions. While 
the ultracold scenario privileges LR interactions, leading to 
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universality in extreme cases|Q], the thermal regime and its 
higher kinetic energies makes (SR) chemical forces prevail. 
As we will show, both SR and LR interactions play the game 
at the concerned kinetic energies. Nevertheless, their roles 
may be independent enough to allow to associate trends in the 
behavior to different regions of the PES, thus offering some 
insight in the underlying dynamics. 

The S^D) + H 2 (X — > SH + H thermal reaction and 
its isotopic variants were the subject of detailed analysis in 
the past [28-39]. Belonging to the class of insertion reac- 
tions D40D, the system is characterized by the presence of a 
deep well (96 kcal mol -1 ) on its singlet ground potential en- 
ergy surface (the l l A' electronic state of H2S) and by a small 
exoergicity (6.90 kcal moL 1 ). Lee and Liu carried out crossed 
molecular beam experiments which were able to conclude the 
insertion mechanism ll3~7l - l39ll . Their measurements motivated 
a number of theoretical studies. Zyubin et al. l28ll . working at 
a multi-reference configuration interaction (MRCI) level with 
multi-configuration self-consistent field (MCSCF) reference 
wave functions, computed the five (1A', 2A', 3A', 1A" and 
2A") potential energy surfaces (PES) which correlate with the 
S('D) + H2 asymptote. The collinear barrier of the first l'A" 
excited PES was found to be 10 kcal moL 1 high. The other 
three surfaces being repulsive, the reaction at moderate ener- 
gies may be assumed to occur on the 1A' PES. An improved 
ground state, based on the same ab initio data set was subse- 
quently produced by Ho et al. B29I1 . This PES is the one we will 
use and it will be labeled as S 1 . Working on it, the quantum 
calculations of Honvault and Launay ID2I1 . precursors to the 
ones we will show below, constituted landmark single-PES 
calculations on atom-diatom insertion reactions at thermal en- 
ergies. Relevant to the discussion of our results will be the 
statistical approaches to the system [33, 34, 41-43]. Let us re- 
call the statistical model of Rackham and Manolopoulos 044ll 
and the quasiclassical statistical model ll45tl . In Refs. ll33ll and 
ll45ll the latter were shown to account very well for the QM 
results by Honvault and Launay. More recently, non-adiabatic 
effects in the product channel were introduced in the quantum 
model in order to calculate reaction cross-sections resolved in 
rotational and fine-structure product states 13511 . 

Much lower in energy in the reagent region there exist elec- 
tronic surfaces which correlate to the S( 3 P), cross the l A' sur- 
face in the H2S well and lead to the same asymptote in the 
product valley. The role of the intersystem crossing in the 
S( 3 P,'D) + H2 collisions was analyzed by Maiti et al. 13611 us- 
ing trajectory surface-hopping methodology. Their calcula- 
tions conclude that the electronic quenching process, S('D) + 
H2 — > S( 3 Po,i,2) + H2, plays a major role in the removal of 
S('D) at energies ~ 250K. The significance of non-adiabatic 
effects in the title system was also revealed in a recent quan- 
tum study performed on the isotopic variant S('D) + HD l46ll . 
We thus expect a significant contribution of the quenching 
process in the current regime. 

In this work we will show the results of new quantum me- 
chanical calculations for the title collision at very low tem- 
peratures. They were obtained starting from an accurate de- 
scription of the LR interactions, essential in order to describe 
collisions at low energy ||47|. In a regime where the kinetic 



energy is so small, LR interactions determine the amount of 
incoming flux which reaches the SR, where rearrangement 
may occur lEoll . In our case the system is characterized by the 
presence of a significant quadrupole-quadrupole (QQ) contri- 
bution, which should lead to important reorientation effects. 
Therefore, the dynamical methodology has been revisited in 
order to include the most of the anisotropy that an adiabatic 
electronic treatment may enable. The theoretical results will 
be compared with the experimental results from the work by 
Costes and co-workers and Sims and co-workers 02711 . 

The paper is structured as follows. In the next section, we 
will briefly describe the new LR calculations, while further 
details will be given in the Appendix. In section III, the hyper- 
spherical approach to the dynamics[48] will be recalled and 
small methodological improvements will be explained. The 
results from the dynamical calculations will be shown in sec- 
tion IV and they will be discussed and compared with the ex- 
perimental data in Section V. Finally, a summary of the work 
and the conclusions will be given in Section VI. 



II. LONG-RANGE INTERACTIONS 
A. Formalism 

The LR interactions between the sulfur atom S('D) and 
hydrogen molecule H2{X 'Sp have been described follow- 
ing perturbation theory up to second order, using a multi- 
polar expansion of the electrostatic interaction operator H e [. 
To account for the 5-fold degeneracy of the open-shell atom 
S( l D), we defined a set of asymptotically degenerate diabatic 
states, taken as a product of the atomic \LA) W and diatomic 
|0) ( - fi - ) unperturbed electronic states. The quantum number 
A = 0, ±l,+2 is the projection of the atomic orbital angu- 
lar momentum L along the Body-Fixed (BF) z-axis, chosen 
along the intermolecular vector R, and the projection A — of 
the diatomic orbital angular momentum relates to the diatom 
axis r, chosen to be in the BF xz-plane. Following previous 
works on open-shell systems [49-52], the matrix elements of 
the first-order and second-order perturbation operators in the 
|L/l) (A) |0) (B) basis set lead to LR potentials that depend on the 
internal coordinates {R,ff), with R being the intermolecular 
separation and the angle between R and r. The dependence 
of the LR potentials on r, the internuclear distance of H2, has 
been neglected in the present case. We have fixed r at its vi- 
brationally averaged value in the ground state, (r),, =0 . 

The anisotropy of the LR S('D)-H2 interaction can induce 
transitions between the fragment states \LA) ( - A) \Q) ( - B) associated 
with different values of the projection quantum number A, 
whereas L is assumed to be conserved throughout the colli- 
sion for large R. We can therefore drop the quantum number 
L and label by \A) the electronic diabatic basis of the S('D)-H2 
interacting system. For linear geometries (0 = 0), there is one- 
to-one correspondence between the states \ A) and the adiabatic 
states of E, n and A symmetry of S('D)-H2. Those adiabatic 
states are labeled according to the value of the quantum num- 
ber A = 0, +1, +2, the projection of the total electronic orbital 
angular momentum _£ along the BF z-axis, with X. being de- 
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fined with respect to the triatomic center-of-mass BF frame. 
For non-linear geometries, the diabatic states \A) are coupled 
through LR interactions, and one has to deal with a 5 x 5 LR 
potential matrix Vaj>(R,0). Since the issue of non-adiabatic 
couplings is not addressed in present work, we chose to keep 
the (complex) signed-/} diabatic basis \A) to express the LR 
potential matrix elements. The diagonalization of the 5x5 
LR potential matrix yields the adiabatic potentials associated 
with the three A' and two A" adiabatic states correlating with 

sOd)-h 2 . 

In the present study, we have considered the electrostatic 
interaction between the permanent quadrupole moments of 
S('D) and FL, leading to electrostatic energies varying as R 5 , 
and the dispersion interaction between the dipole-induced mo- 
ments of the two species, giving rise to dispersion energies 
varying as R . The contribution of induction effects to the 
LR potentials has been neglected, since the dominant in- 
duction contribution involves interactions between permanent 
quadrupole moments and dipole-induced moments, which are 
proportional to R~°. According to Eqs. [TBI and [TBI of the Ap- 
pendix, the LR potential matrix elements write: 

Vw(R,B) = 4 C f C 2 ,a>-a(0,0)--^ J] C K a--a(0,0) 

£=0.2 

(1) 

where the angular functions C/ >m (#, <p) are normalized spher- 
ical harmonics lT76ll . and the interaction coefficients Ci A ' and 
C AA k stands for the electrostatic (quadrupole-quadrupole) and 
dispersion (dipole-induced dipole-induced) contributions, re- 
spectively. 

B. Quadrupole moments and dipole polarizabilities 

The evaluation of the LR electrostatic C AA and disper- 
sion C AA k coefficients requires to determine the permanent 
quadrupole moments and dynamic dipole polarizabilities of 
S('D) and H2. For H2, we have selected accurate values 
from the literature, corresponding with a H2 geometry at 
its vibrationally averaged value, (r),, =0 = 1.449 «o- We 
used the permanent quadrupole value (0|Q2o|0> = +0.481 
a.u. of Ref. 15311 . obtained from full Configuration In- 
teraction (CI) calculations with a large Gaussian-type or- 
bitals basis set, including bond-centered polarization func- 
tions. For the dynamic dipole polarizabilities °°a zz (ia)) and 
00 0^(20/), we chose the values tabulated in Ref. i54ll . obtained 
from a sum-over-states formalism with explicitely electron- 
correlated wavefunctions. The corresponding static polariz- 
abilities, °V,(0) = 6.721 a.u. and °°a xx (0) = 4.739 a.u., are 
in good agreement with other literature values ll53ll . 

For S( l D), we performed ab initio calculations by means 
of the Dalton quantum chemistry code [55]. The permanent 
quadrupole moments have been calculated as the expectation 
values of cartesian quadrupole moment operators, ©„„ with 
u - {x,y,z}, using ab initio MCSCF electronic wavefunc- 
tions. The MCSCF wavefunctions were generated by dis- 
tributing six electrons among 13 orbitals (3 s, 3p, 3d, 4s and 



Ap), the inner-shell orbitals being kept frozen in their ROHF 
form, and the Sadlej-pVTZ basis set was used. The result- 
ing value of the quadrupole moment (LO|22o|W) = +2.075 
a.u. was employed to get the whole set of nonzero matrix 
elements from the Wigner-Eckart theorem. Notice that the 
latter value is identical to the one already reported in the lit- 
erature ll56l for the same level of calculation. The dy- 
namic dipole polarizabilities of S( l D) have been determined 
by means of MCSCF linear response calculations, using the 
MCSCF wavefunctions previously defined. For each substate 
of definite symmetry, we obtained a set of Cauchy moments, 
each of them corresponding with a given cartesian compo- 
nent of the electric dipole operator fi„=x,y,z- We then employed 
analytical continuation techniques following the [n,n- l] a 
and [n,n — 1]^ Pade approximants procedure of Ref. 15811 to 
get lower and upper bounds to the dynamic polarizabilities 
n 'a m (iu>), where v stands for an atomic substate of definite 
symmetry. For n = 9, the dispersion coefficients associated to 
lower and upper bounds of the polarizabilities are converged 
within less than 0.02%. The values of the corresponding static 
polarizabilities are identical to those already published lf56ll . 
From the diagonal elements of the dynamic polarizabilities 
n 'a m (ia>), we derived the whole set of coupled spherical dy- 
namic polarizabilities AA a^i)kq following the procedure de- 
scribed in Ref . JH . 

The resulting LR interaction coefficients for the electro- 
static C AA and dispersion C'l \ contributions are tabulated in 
the electronic supplementary information (ESI). To the best 
of our knowledge, no LR coefficients had been determined 
so far for the S( 1 D)-H2 system, and thus it is not possible 
to compare with the results of previous works. By sum- 
ming the tabulated values over the diagonal matrix elements 
(A = A') we can retrieve state-averaged interaction coeffi- 
cients. The state-averaged electrostatic contribution vanishes, 
and we get for the isotropic and anisotropic dispersion contri- 
butions C(,fl - 40.338 a.u. and C6,2 = 4.146 a.u., respectively. 



III. DYNAMICAL METHODOLOGY 

A. The hypersperical approach 

The quantum methodology used to carry out the dynami- 
cal calculations was described in previous works on alkalis at 
ultracold energies 05 9l 16011 . and elsewhere in the context of 
thermal reactive scattering[61]. In fact, most of the conver- 
gence parameters used in this work are the same which were 
used in the study of the title collision at thermal energies ll3lll . 
Let us simply recall that in the hyperspherical quantum re- 
active scattering method developed by J. M. Launay|48J the 
configuration space is divided into inner and outer regions. 
The positions of the nuclei in the inner region are described in 
terms of hyperspherical democratic coordinates. The logarith- 
mic derivative of the wavefunction is propagated outwards on 
a single adiabatic PES. At a large enough value of the hyper- 
radius the former is matched to a set of suitable functions, 
called asymptotic functions, to yield the scattering S-matrix. 
We chose an intermolecular separation of * 10 a.u. for the 
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matching. 

The asymptotic functions provide the collisional boundary 
conditions. When working at thermal energies they are fa- 
miliar regular and irregular radial Bessel functions which ac- 
count for the presence of the centrifugal potential at large in- 
termolecular separation. Recently, they were propertly mod- 
ified to include the effect of an isotropic dispersion interac- 
tion (~ R~ 6 ), thus enabling the study of ultracold collisions 
of alkalis lf59l l60h . In the current work further changes are re- 
quired. Asymptotic functions must account also for the pres- 
ence of the anisotropic electrostatic and dispersion LR inter- 
action, introducing reorientation effects which may take place 
while the reagents approach. More details about their defini- 
tion are given in the following section. 



B. Introducing the anisotropy in the outer region: the 
asymptotic functions 

Using the set of Jacobi coordinates [R, r, 9) corresponding 
to the S+H2 arrangement, the total Hamiltonian of the system 
can be expressed as 



H 2 1 d 2 

H = =R + 

2// R dR 2 



1 



n 2 1 d 2 



1 



Y-— - — r + 

2jxR 2 2m r dr 2 2mr 2 



f + H d (2) 



where H e i is the electronic Hamiltonian, 1 is the orbital angu- 
lar momentum of the atom with respect to the center of mass 
of the diatom, and j is the rotational angular momentum of the 
latter. The total angular momentum of the nuclei (conserved 
in an adiabatic approach) is given by J = j + 1. In the current 
study, we chose an adiabatic treatment of the dynamics, as- 
suming that the collision occurs only on the ground adiabatic 
PES. Hence, among the five adiabatic PESs correlating with 
S('D)+H2 (three A' and two A" singlet states), we consider 
only the lowest one, associated with the ground 1 l A' elec- 
tronic state. Hereafter, the lowest eigenvalue of H e i will be 
labeled by V°(R, r, 9). The diagonalization of the LR poten- 
tial matrix of Eq. Q]provides a value for such energy, accurate 
at large intermolecular separations, with r fixed to its vibra- 
tional^ averaged value (r) v=0 . 

A convenient basis in order to expand the nuclear wave- 
function in the LR region is the one characterized by quantum 
numbers (J, M, v, j, /), with (v, j) the rovibrational quantum 
numbers of the diatom, / the relative orbital angular momen- 
tum and (J, M) the total angular momentum and its projection 
on the Space-Fixed Z axis. We will call this basis Bl, and rep- 
resent it as (f>™. Such a basis, well adapted to handle Coriolis 
couplings, is used in the hyperspherical approach to expand 
the asymptotic wavefunctions, which are matched with the SR 
information obtained in hyperspherical coordinates. 

If the system approaches collision with quantum numbers 
(J, M, vo, jo, h), we will asumme that (in addition to J and M) 
the rovibrational quantum numbers, (vo, jo), remain well con- 
served in the LR region. This is well justified in the present 
case given the large energy gap between different rovibra- 
tional states relative to the considered collision energy range. 
Within this approximation the total electro-nuclear wavefun- 
cion, X ¥ JM , , can be expanded in the LR region as 
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F, (R) 



R 



(3) 



where |1 l A') is the ground adiabatic electronic state and all 
the "conserved" quantum numbers have been suppressed in 
the notation of the radial coefficients, f'°(R). Introducing the 
expansion (0) into the time-independent Schrodinger equation 
associated with a total energy E, EN 1 = and using the 
Hamiltonian in Eq. |2j it is straightforward to obtain the fol- 
lowing system of coupled radial equations 



n 2 d 2 1(1 + l)tr 
'^dR 2 + 2fiR 2 



F' l °(R) = -J](V) u ,(R)F ! ?(R)(4) 



where E c , the collision energy, is given by E - E v j, with E v j 
the internal energy of the diatom. (V)ij>(R) designs the matrix 
elements in basis Bl of V°(R,r,ff) - V Hl (r), with V# 2 (r) the 
asymptotic H2 diatomic potential. By inwards integration of 
Eq.0we obtain the "regular" (Ff )k ) and "irregular" (F^' ) 
asymptotic radial wavefunctions, corresponding to an incom- 
ing (J, M, vo, 70, '0) channel. They are defined as the ones 
which behave as 



fi->co 
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sm(kR-ln/2)6i h /k 
■cos(kR-lnll)6 lh lk 



1/2 



(5) 
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where k is the wavenumber of the considered channel. Such 
radial functions account for the correct LR behaviour to match 
with the inner-dynamics in the presence of an anisotropic 
V°(R, r, ff) potential. The coupled-equation version of the 
method of De Vogelaere [62] is used to solve for them. 

Regarding the calculation of the potential matrix elements 
(V)ij'(R), it is convenient to define an intermediate nuclear 
basis, B2, labelled by the projection Q.j of j and J on the BF 
z-axis. The basis set B2 is given by 



,JM 



^^0&/ar,Ay)r A (0,O), (6) 



where D J ^ n designs a Wigner rotation matrix element and 
(a,/3, y) are the Euler angles corresponding to the transforma- 
tion between SF and BF frames. The change of basis between 
Bl and B2 simply involves 3j symbols. The matrix elements 
of the ground PES in basis B2 are given by 



<V> ; 



ijSi>(R) = 5n,n'2n J 4/r)I%.(0, 0)V°(R, r, 9) sin 9 dr d9 

(7) 

The potential does not couple states with different Qy, al- 
though reorientation effects in the BF system happen due to 
the Coriolis term. Assuming in Eq. |7]that the variations with 
R of the mean value of V°(R, r, 9) in the rovibrational state 
Xv,ff) are we H approximated by the variations of the potential 
at the averaged bond distance, r = (r) v j (assumption which 
has been proven valid for the current PES), we reach the con- 
venient expression 
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(8) 



Once the potential matrix elements are calculated in basis B2 
we change to basis Bl, thus obtaining the elements (V)ij>(R) 
involved in Eq. @] 
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<V>Q^(*). 



Let us finally note that {V)i oJo (R) + k(k + l)h 2 /2/iR 2 has the 
meaning of the effective potential felt by the colliding partners 
at a distance R when approaching in the state <f> J v ^ i - Such 
meaning will be used below. 



C. The potential energy surface 

In our adiabatic approach we consider that the collision 
takes place on the ground PES. The global fit for the l l A' 
PES performed by Ho et al. H29fl (surface SI) was widely used 
in the past and its SR region was tested by comparison with 
experiments at higher energies I30ll3lll . Nevertheless, the sur- 
face does not pay the necessary attention to the description of 
the LR region, that being the reason why we performed our 
own calculations of the LR potentials (see Sec. [III). To com- 
bine the SR interactions given by S 1 with our description of 
the LR region (the lowest eigenvalue of the LR potential ma- 
trix in Eq. [TJ, we performed a smooth switching from one 
to the other. Thanks to the approximation in Eq. [8] only the 
LR potentials for fixed r = (r) v j ~ (r) v =o are required, and 
only the switching for such internuclear distance is consid- 
ered. The switching was performed around the central value 
Ro = 12.5 a.u., within a range of intermolecular distances 
given by AR = 3 a.u. The angle-independent switching func- 
tion f(R) = {cos[tt(R + AR/2 - R )] + l}/2 was used. In con- 
trast, note that an intermolecular separation of » 10 a.u. was 
considered for the matching of the inner dynamics with the 
external asymptotic functions. 

It is worth mentioning that a new global 1 l A' PES for the ti- 
tle system was recently published B63I1 . A term to account for 
the (~ R~ 6 ) dipole-induced dipole-induced (DD) dispersion 
interaction was included in the global fit. Nevertheless, the 
quadrupole-quadrupole (QQ) electrostatic interaction (~ R ), 
which dominates the LR behavior, was not explicitly consid- 
ered. As such term is relevant to describe the reactivity of 
orf/io-hydrogen at low collision energies, we have preferred 
to employ the PES obtained from the switching procedure just 
described. 



D. The collision with j =0 para-hydrogen 

The collision of S('D) with para-hydrogen in its ground 
rotational state deserves special attention. Within a 
non-adiabatic treatment the effective contribution of the 
anisotropic terms in V^viR, ff) (Eq.Q} is found to vanish. This 
is not the case when using an adiabatic approach instead. In 
order to show this we define a set of electro-nuclear basis 
functions in terms of the diabatic electronic basis defined in 
SectionHIl 



,J T M T 



Xv,j(r) 



^LllD^^fi, y)Y j a ] (6, 0)1^X10) 



where Jj designs a total angular momentum, Jt = j + 1 + L, 
which in contrast with Eq. [6] includes the atomic electronic 
angular momentum, and Q7- (= Clj + A) and M T label its pro- 
jection on the BF and SF z-axis, respectively. The use of this 
basis in combination with the whole interaction matrix given 
by Eq. [TJ would efficiently account for the exchange of an- 
gular momentum among electrons and nuclei (Jt being con- 
served) and reorientation effects appearing while the reagents 
approach ll64ll . Potential couplings between two different 
channels in this non-adiabatic treatment would involve the in- 
tegration in 6 of the product of two spherical harmonics (orig- 
inating from Eq. [Tot and one matrix element Vaj?(R,0). The 
latter consists for spherical harmonics, Yf, m , with k = 0, 1,2 
(Ck, m ~ Yk, m in Eq. [TJ. When considering the case of para- 
hydrogen with j=0, such integral is found to vanish for k + 0. 
The contribution of QQ and dispersion anisotropic terms dis- 
appears, and thus the R 5 dependence. Only the isotropic dis- 
persion term, ~ R~ 6 , is found to play a role. This is not the 
case when the adiabatic approach is used. The diagonaliza- 
tion of the LR V,\j>(R,ff) matrix yields a ground PES with a 
R~ 5 dependence, which is not cancelled through the integral 
in Eq. [8] In order to correct such a drawback of the adiabatic 
treatment we artificially matched S 1 with a pure isotropic dis- 
persion term, -C(,/R 6 , when dealing with j=0 para-hydrogen. 
In this way we introduced in the adiabatic treatment the effec- 
tive LR interaction which the colliding partners should expe- 
rience. The C 6 '? coefficient, the lowest isotropic one, was 
chosen for that purpose. 



IV. RESULTS 

Quantum reactive cross-sections for collisions of S('D) 
with j=0 para-hydrogen (p-Fb) and j=l orf/10-hydrogen (o- 
H2), in the energy range 1-120 K, were calculated using the 
methodology described in the previous section. In Fig. [JJ the 
obtained results are shown. A small dependence with the ro- 
tational state and some interesting oscillations, more apparent 
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FIG. 1: Quantum reaction cross-sections for the collision of S('D) 
with j=0 P-H2 and j = l 0-H2 at low kinetic energies. 



in the case of p-H2, are relevant features. The cross-sections 
were averaged using the room-temperature ratio 3:1 of 0-H2 
to /7-H2, so-called normaZ-hydrogen (71-H2), which is expected 
to be valid in the conditions of the experiment by Costes and 
co-workers 

IH. They were also convoluted with a Gaussian 
distribution of collision energies to account for the velocity 
and crossing angle spread in the experiment [27]. In Fig. [2 
the experimental excitation functions (in arbitrary units) were 
scaled to minimize the sum of the squares of the differences 
between every experimental result and its corresponding the- 
oretical value. Both sets show very similar decreasing trends 
as a function of the collision energy. 

Some additional calculations for particular collision ener- 
gies in the range 150- 1 500K were performed in order to allow 
interpolations and the calculation of thermally-averaged rate 
constants to compare with the experimental results by Sims 
and co-workers 1231 12711 ■ Again, the room-temperature ratio 
3:1 was considered. It was estimated that such proportion is 
not perturbed by the rapid cooling to 77 K and the subsequent 
adiabatic expansion to 5.8 K in the experiment. By using a 
Maxwell-Boltzmann distribution the population of higher ro- 
tational states was estimated to be negligible for all the ex- 
perimental temperatures except for the highest one, T=300K. 
Given the independency of the sets of rotational states with 
even and odd quantum number and the big rotational constant 
of H2 (~ 85K), only in the case of p-Bfe and for such high 
temperature the population of j = 2 was found to be sim- 
ilar to the population of the ground state. Accordingly, we 
calculated the cross-sections for j = 2 at some particular col- 
lision energies and concluded a difference of ~ 5% between 
the cross-sections for j ' — and j = 2. The contribution of 
j = 2 was not considered in the calculation of thermally aver- 
aged rate-constants and a value somewhat lower than the one 
in Fig.[3]is expected for T=300K. 

Interestingly, the theoretical rate constants approximate 
fairly well the experimental ones, remaining below them with 
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FIG. 2: Comparison of the excitation function obtained in the exper- 
iments by Costes and co-workers [27] with our theoretical results. 
The latter were averaged using the ratio 3:1 of 0-H2 to p-Yio, and 
convoluted with a Gaussian distribution of collision energies to ac- 
count for the velocity and angle spread in the experiment |27|]. The 
excitation function, in arbitrary units, was scaled in order to optimize 
the agreement with the theoretical data (see text). (Experimental data 
reproduced with the permission of the authors). 



the exception of the value for the lowest experimental tem- 
perature (5.8K). This temperature is also the most challenging 
from an experimental point of view. The 300 K experimental 
result of Black and Jusinski[65] is also shown for compari- 
son, in perfect agreement with the recent measurements. Let 
us remark that a multiplicative factor 1/5 was included in the 
calculation of rate constants, as only one of the five asymptot- 
ically degenerate adiabatic surfaces is assumed to lead to the 
products in the considered energy range. 



V. DISCUSSION 
A. The significance of the long-range interactions 

In order to stress the need for accuracy in the description 
of LR effects while working at low collision energies, the 
sensitivity to such effects is shown in Fig. [4] Cross-sections 
for the collision with /7-H2, calculated using an uncorrected 
PES SI [29], are compared to the ones obtained by matching 
it with a pure isotropic ~ R 6 dispersion term at LR (see Sec- 
tion UnDji. As the figure shows, for energies below 20K the 
accurate inclusion of the LR interactions appears as essential. 
The use of S 1 without corrections leads to differences of the 
order of 30% for energies around IK. 

Note that the LR behaviour of an accurate ab initio surface 
for the ground electronic state contains a non-negligible R 5 
character in the title system. As argued above, such contri- 
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FIG. 3: Comparison of the experimental total removal rates (filled 
circles), obtained in the experiments by Sims and co-workers 123L 12711 
with the thermally averaged reaction rate constant obtained in our 
calculations (dashed line). Again a ratio 3:1 of 0-H2 to P-H2 was 
assumed. The 300 K experimental result of Black and Jusinski[65] 
(marked with a cross) is also shown. (Experimental data reproduced 
with the permission of the authors). 
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FIG. 4: Sensitivity of the calculated cross-sections to the LR intera- 
tions. Cross-sections for the collision of S('D) with p-H 2 , obtained 
by using an uncorrected PES SI |29] are compared to the ones ob- 
tained by matching it with a pure isotropic ~ R~ 6 dispersion term, 
thus accounting for the expected LR effective behaviour. 



bution would effectively cancel while using a non-adiabatic 
approach (which considered the 5 asymptotically degenerate 
surfaces and all the non-adiabatic couplings). Nevertheless, 
it does not cancel while using an adiabatic approach (that be- 
ing the reason to force a LR behaviour R in the calculations 
with /7-H2). The need for such artificial modifications rein- 
forces the idea that for asymptotically degenerate systems an 
accurate treatment of the LR non-adiabatic dynamics seems 



essential. To the best of our knowledge, non-adiabatic effects 
in low energy collisions have not been explored in depth up to 
the date. 



B. The electronic quenching 

In view of the results of Maiti et a/.j36| and Tian-Shu et 
a/.|H3l the electronic quenching is expected to play a role in 
the title collision. Since the experimental rate constants mea- 
sured by Sims and co-workers 112711 probe the total removal 
of S( D) by n-H2, they should comprise the disappearance of 
S( l D) due to the reactive process (formation of SH+H) or to 
the electronic quenching (relaxation to S( 3 Po,i.2))- The ex- 
citation functions measured by Costes and co-workers ll27ll 
monitor only the first mechanism. As our theoretical calcula- 
tions do not include the possibility of quenching some com- 
ments are pertinent in order to justify mutual comparisons. 
What should we expect from the use of an electronically- 
uncoupled (EU) formalism to describe the dynamics when an 
electronically-coupled (EC) one would be more appropriate? 
A reasoning according to the statistical model of Rackham 
and Manolopoulos[44], applied in the past to the system at 
thermal energies 113 3i 14211 . may help to clarify the limits and 
validity of our approach. If the collisional process is assumed 
to be mediated by a long-lived complex (due to the deep well 
in the PES) we can distinguish the step of formation from the 
step of statistical decomposition of the complexes. Since the 
singlet-triplet crossing occurs inside the well we will consider 
the quenching as an alternative way of decomposition. Be- 
sides, we affirm that our calculation describes correctly the 
first step, the way the probability is captured. Regarding the 
second step, in a scenario where some complexes would de- 
compose to give back the reactants (I), some to give reac- 
tion products (R), and the rest would suffer quenching (Q), 
our electronically-uncoupled (EU) calculation does not allow 
this latest possibility. By making use of the principles of the 
model, and designing by crf c , cr^ c , <x^ c and cr^ u , cr^ u the 
cross-sections for each possible outcome in an hypothetical 
EC (unprimmed) or in our EU (primmed) calculation, we will 
show that 

<4 c + cjf><4 u >oF (id 

As a consequence, the theoretical rate constant calculated us- 
ing <t^ u constitutes a lower bound for the experimental to- 
tal removal rate of S{ l D). The latter should be given by the 
sum of both reactive and quenching contributions in an exact 
EC calculation. Besides, cr^ u is itself an upper bound for the 
real reaction cross-section, monitored in arbitrary units by the 
measured excitation functions. These statements should be 
valid even if <r% u is far from cr^ c . Below we will show how 
to reach these conclusions. 

According to the statistical model l44ll . the fraction of 
complexes which decompose into a particular channel is re- 
lated to the capture probability of forming the complex start- 
ing from that particular channel. Let us denote with A J n (E), 
B J n ,(E) and C J n „{E) the energy dependent capture probabili- 
ties starting from reagent, product and quenching sides, cor- 
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responding to particular channels n,n',n", respectively. J 
labels the total angular momentum. Let us define the sum 
in open channels of the capture probabilities associated to 
each possible outcome, A J (E) = £ A J n (E), B J (E) = £ B J n ,(E) 
and C J (E) = Y*C J n „{E). Following the prescriptions of the 
model, the reaction probability corresponding to the reagents 
colliding in channel n, P Rn (E), is given by P R EU {E) = 
A ] n {E)B ] {E)l(A J (E) + B J (E)) in our EU calculation and by 
P R EC (E) = A J n (E)B J (E)/(A J (E) + B\E) + C\E)) in an hy- 
pothetical EC calculation. Interestingly, the former is larger. 
Now, essentially by summing over J, one obtains cross- 
sections and it follows that o- R u > cr R c , which is the sec- 
ond inequality in expression [TT] A similar reasoning leads 
from the probability of total removal, P J R EC {E) + Pg E n c (E) = 
A J n (E)(B J (E) + C J (E))/(A J (E) + B J (E) + C J (E)\ to the first 
inequality. In agreement with it, the thermally averaged the- 
oretical rates in Fig. [3] remain below the experimental ones, 
with the exception of the value for the lowest temperature. 



C. The Langevin-type capture dynamics 

Both the theoretical and the experimental excitation func- 
tions, in Fig. |2] show no threshold, consequence of the attrac- 
tive character without barrier of the reaction on the ground 
electronic state. Remarkably, both follow very similar trends 
as a function of the kinetic energy, in spite of the absence of 
quenching channels in the theoretical calculations. As a first 
guess, such proportionality would result from the real reaction 
probabilities following the entrance channel capture probabil- 
ities, A J n {E), given that the capture is well described in our ap- 
proximate calculation. Moreover, if our guess were correct the 
experimental reaction cross-sections should have a Langevin- 
capture character l66ll (see note ll77lo . 

Such qualitative reasoning can be put in solid grounds by 
making use again of the statistical model. Given the exother- 
micity of both the reactive and quenching processes, many 
internal states are open even at low collision energies in the 
corresponding arrangements. In particular, B J (E) may be 
approximated by the number of open channels, thus being 
much larger than A J (E). The value of C J (E) is more diffi- 
cult to estimate as the capture from the quenching side re- 
quires a non-adiabatic transition, of unknown probability, to 
occur. Regardless, it is resonable to accept that B J (E), and 
thus B J (E) + C J (E), are much bigger than A J (E) in the en- 
ergy and angular momentum ranges which contribute, thus 
following that most of the complexes are expected to form 
SH+H or to suffer quenching. Besides, such figures are ex- 
pected to be essentially constant, B J (E) ~ B, B J (E) + C J (E) ~ 
B + C. Accordingly, in we neglect the contribution of A J (E) 
in the denominator of P R EU (E) (see above) we find a re- 
sult which is essentially equal to a reactant capture proba- 
bility, A J n (E). Proceeding similarly we find that P R EC (E) 
is a constant portion of the capture probability, P R EC (E) ~ 
(B/(B + C))A ] „(E), and thus proportional to P J R EU (E). Finally, 
the sum P R EC (E) + Pg EC (E) is again equal to the capture 
probability A J n (E) and therefore to P R EU (E). By summing 



over J to obtain cross-sections, the previous results are found 
equivalent to the following statements: i) cr R u is essentially a 
capture cross-section which accounts for the total number of 
complexes formed in the collision, ii) cr R u should be found 
proportional to cr R c , and iii) o- EC + ct E q and o- R u must be 
similar to each other. We will check their validity below. 

Regarding the first statement, the quantum cross-section 
averaged for n-Wz is compared in Fig. [5] with an averaged 
Langevin capture estimate and a good agreement is found. 
A capture dispersion potential -C^IR 6 was assumed for p- 
H2, what leads to an analytical Langevin cross-section given 
by o-^iE) = 37r(C 6 /4£ , ) I/3 . A QQ electrostatic poten- 
tial, -C5/R 5 , and the corresponding cross-section, o- 5 Lan (E) = 

(C 5 /6V3£) , were assumed for 0-H2 instead. The val- 
ues C 6 = C°'£ =0 = 38.9 a.u. and C 5 = 4.08 <e20>s«220>ff 2 
were respectively taken. The former is the lowest isotropic 
dispersion coefficient and the latter results from an angular 
average of the ground eigenvalue of the QQ contribution in 
Eq. Q] Double weight is assigned to the perpendicular ap- 
proach, C 5 (6 = 90) = V39/2 <e20>s<{22o>H 2 , with respect 
to the linear one, C$(6 = 0) = 6 (620)5 <fi2o)// 2 . the fac- 
tor 4.08 resulting from such weighted average. The second 
statement explains the similar trends found between theoreti- 
cal and experimental excitation functions, in spite of the non- 
inclusion of quenching in the calculations. Finally, the third 
statement justifies the similarity of the experimental and the- 
oretical rates, shown in Fig. [3] The bigger discrepancy be- 
tween the experimental and theoretical rate coefficients at 5.8 
K, the latter being larger, may be related with the neglected 
non-adiabatic couplings between the three singlet A' states 
which correlate with the reactants. In fact, the discrepancy can 
be interpreted as the theoretical simulations capturing more 
probability in the complex than they should, and neglecting 
the couplings to repulsive states could be a reason for such an 
effect. The effect of such couplings should increase dramati- 
cally when decreasing the kinetic energy, while the contribu- 
tion from the neglected singlet-triplet couplings is expected to 
be less energy dependent. Regardless, such low kinetic ener- 
gies should be more sensitive to possible inaccuracies of the 
considered PES. 

To outline this section, let us remark that our theoretical 
cross-sections are essentially capture cross-sections and they 
must account well for the total number of complexes formed 
in the collision. The number of complexes which eventu- 
ally react is found a constant portion of the total number. 
This explains the experimental reaction cross-sections being 
proportional to the theoretical cross-sections. Finally, as the 
experimental total removal rates include the sum of reaction 
and quenching, they are monitoring again the total number of 
complexes, so they are well approximated by our theoretical 
results. 



D. The structures in the cross-sections 

The theoretical cross-sections show interesting oscillations 
superimposed to the overall Langevin behaviour (Fig. |5J, 
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FIG. 5: In the figure, the theoretical cross-section averaged using 
the 3:1 proportion of o-H 2 to p-H 2 present in the experiment by 
Costes and co-workers 1 27] is compared with the averaged analyt- 
ical Langevin estimate|66] (see text). 



They can be observed for both p-H2 and 0-H2, although their 
presence seems to be washed out by the sum in partial waves 
in the latter case. In a naive approach, zero-order resonances 
present in the case of p-H2 would split in the case of 0-H2 
due to the different values of Q. Their abundance and mutual 
overlap would justify their being more difficult to observe in 
the latter case. 

The observed structures, which can be attributed to partic- 
ular partial waves and channels have a double origin: some of 
them display the effective opening of a partial-wave /; others 
seem to be due to the presence of shape resonances. Regard- 
ing the former, the Wigner threshold laws provide a behaviour 
for the I cross-section below the centrifugal barrier given by 
£/-i/2 jf^ a f ter j-jjg Carrier were overcome, the reaction prob- 
ability reached quickly the limit value of one, the behaviour 
should follow a 1/k 2 dependence, that is ~ £'~ 1 ll67l469ll . This 
would result as a "bump" in the cross-sections. A bump of 
this type should happen at energies over the height of the cen- 
trifugal barrier. On the contrary, oscillations for high partial 
waves steming from shape resonances would occur at energies 
below the centrifugal barriers, the associated effective poten- 
tials providing some kind of trapping. 

Some structures have been selected and they are shown in 
Fig. [6] to be associated to particular partial waves and chan- 
nels, <p\!fi- The corresponding centrifugal barriers are also dis- 
played as vertical lines. They were calculated using effec- 
tive potentials, (V)u(R) + 1(1 + l)h 2 /2/iR 2 , where the average 
(V)ij(R) was defined in Eq. [9] Panel (a) corresponds to j=0 
P-H2 and panel (b) to j=l 0-H2. Given the shape and location 
(above the centrifugal barrier) of the peaks associated to /7-H2, 
with entrance channels J — I = 7 and J - I - 12, and to 0-H2 
with entrance channel J — 9, 1 = 8, we conclude that their 
origin is the opening of a new partial wave. In contrast, the 
peaks labelled A and C for P-H2, associated to the J - I - 15 




Collision energy (K) 

FIG. 6: In the figure, some oscillation observed in the quantum re- 
action cross-sections are attributed to particular partial waves and 
entrance channels. Panel (a) corresponds to j=0 p-Ha and panel (b) 
to j=l 0-H2. Peaks associated to J = I = 7 and J = I = 12 in the case 
of P-H2 , and J = 9, 1 = 8 in the case of 0-H2 occur at energies just 
above the corresponding centrifugal barriers (shown as vertical lines 
in the same color or type of line). In contrast, peaks A, B, C and D 
occur below barrier (see text). 



and J — I — 16 contributions, are found below barrier. The 
same happens for the peaks A, B, C and D in the lower pan- 
nel, associated to the entrance in channels [J = 16, 1 = 15), 
(J = 16,1= 16), (J = 17,1 = 16) and (J = 11,1= 17), re- 
spectively. As we will explain below, we attribute them to the 
presence of shape resonances. Note that peaks A and C are 
found at similar energies for 0-H2 and p-fy. 

The presence of multiple resonances in the reaction proba- 
bility of the title system for zero total angular momentum was 
reported in the past 1 30] . They manifest the indirect mecha- 
nism through the formation of long-lived complexes. If the 
peaks selected in the current work were resonances, the as- 
sociated trapping should happen in an intermediate range of 
radial distances. On the one hand, they appeared in calcula- 
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FIG. 7: In the figure, the effective potentials, (V),j(R) + 1(1 + 
VfiP /2/uR 2 , of the channels associated to peaks A, B, C and D (see 
Fig. [6} are shown to provide the necessary trapping to support quasi- 
bound states in the form of shape resonances. Panel (a) corresponds 
to j=0 p-H.2 and panel (b) to j= 1 0-H2. The horizontal lines mark the 
location in energy of the peaks. Note the similarity of the effective 
potentials corresponding to the pair of channels (J = 16, / = 15), 
(/ = 17, / = 16) of o-H 2 to the ones corresponding to J = I = 15, 
J = I = 16 ofp-H 2 . . 



tions where the LR interaction was not well accounted for (see 
Fig. |4j, so they cannot be located at very large distances. On 
the other, their location at very short range would hamper, due 
to the high anisotropy of SR chemical forces, any rotational 
adiabaticity which may lead to similar effective potentials, 
and thus to the similar structures observed, for j=0 and j=l. 
In good agreement with such qualitative reasoning, the PES 
shows a shallow well in the entrance valley at distances of ~ 7 
a.u. ll42ll - With a depth of ~ 500K for collinear approaches, 
and appearing as a small secondary minimum opening the ac- 
cess to the big well for insertion paths, this complex is the 
probable origin of many of the structures in the cross-sections. 
As a result of such well, the effective potentials of the channels 



associated to the peaks show a double maximum structure (see 
Fig-El- The local minimum inbetween may support shape res- 
onances at precisely the energies where the structures appear. 
Remarkably, the effective potentials corresponding to chan- 
nels (J = 16,1 = 15) and [J = 17, 1 = 16) of 0-H2 are respec- 
tively very similar to the ones corresponding to J = / = 15 and 
J = / = 16 of p-H2- This justifies the appearance of peaks A 
and C at similar energies in both cases. The approaching part- 
ners would be trapped in this preliminar well until they tun- 
neled into the big one, thus finding the way to reaction, or they 
tunneled backwards, thus giving the reagents back. The qua- 
sibound states seem to prefer "prereaction"though, increasing 
reactivity at the associated energies. Note finally that Fano- 
type profiles, which indicate quick changes in the phase of 
the elastic S-matrix elements, are found in the elastic cross- 
sections at similar energies, thus reinforcing the hypothesis 
about the nature of the peaks as resonances. 

Let us remark that the role of van der Waals wells in reactive 
collisions at very low energy was recently stressed B70I1 . As 
stated, the supported quasibound states may lead to a tempo- 
ral trapping of the system at short distances, and to an increase 
of tunneling through the barrier. Such effects were mainly dis- 
cussed for tunneling-dominated reactions like Li+HF H7 111 and 
Cl+HD [70], where the reaction would be closed otherwise. 
In addition to the considered shape resonances, Feshbach-type 
ones may occur for zero angular momentum and can dramat- 
ically modifying ultracold reaction cross-sections and zero- 
temperature rate constants ll72tl . 



VI. CONCLUSIONS 

In this work, the S('D) + H2 collision was analyzed in 
the kinetic energy range 1-120K. Within an adiabatic treat- 
ment, working on the lowest surface which correlates with the 
reagents, cross-sections and thermally averaged rate constant 
were obtained using the quantum reactive scattering method 
developed by Launay et al. 14811 . The resulting picture of the 
process is the following: the behavior can be understood on 
the average in terms of a classical Langevin (capture) model 
where centrifugal barriers determine the amount of reactive 
flux which reaches the barrierless transition state; addition- 
ally, the structure of the van der Waals well provides tem- 
poral trapping at short distances, thus helping the system to 
find its way to reaction for high partial waves. Comparison 
with landmark experimental data obtained using the crossed 
molecular beam machine with variable beam intersection- 
angle JH (Bordeaux) and the CRESU technique JH |H 
(Rennes) shows a fairly good agreement. Accurate electronic 
calculations of the LR interactions were required in order to 
propertly describe collisions at such low kinetic energy. In 
particular, the anisotropic QQ term of the interaction appears 
to be important to describe the collision with 0-H2, although 
it does not contribute in the case of p-FL;. 

Previous theoretical works at higher energies indicated that 
the electronic quenching process, and thus the intersystem 
crossi ng, may play a significant role in the total removal of 
S('D)|36, 46]. Nevertheless, a good agreement of our theoret- 
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ical and the experimental total removal rates has been found 
in spite of the adiabatic treatment of the collision. More- 
over, the reescaled experimental excitation functions follow 
very well the theoretical ones. An small contribution of the 
quenching must not be concluded though. The current the- 
ory accounts well for the total number of complexes formed 
in the collision, this amount being what our theoretical rates 
and cross-sections essentially monitor. Now, as most of the 
complexes decompose to give either reaction or quenching, 
and experimental removal rates sum up both contributions, 
a good agreement of theoretical and removal rates follows. 
Besides, the ratio of reaction to quenching is not expected 
to significantly change in the small considered energy range. 
The amount of complexes which 'react' at one particular en- 
ergy appears then a constant fraction of the total number of 
complexes captured at such energy. A proportionality of the 
'real' reaction cross-section to our theoretical cross-sections 
automatically follows, and therefore the good agreement with 
the experimental excitation functions. It must be noted that 
both discussed agreements seem to follow independently to 
the proportion of quenching and, in our opinion, the question 
remains open. Future measurements of the branching ratio 
reaction/quenching would be thus desirable. 

Interestingly, the theoretical cross-sections showed oscillat- 
ing features. We conclude that some of them are associated to 
the presence of a small van der Waals well in the entrance val- 
ley of the ab initio PES SI 12911 . Far from being an artifact 
of the fit, this well seems to be real, appearing in some ab 
initio calculations we have performed to verify its existence. 
Nevertheless, we have found a much smaller depth. Although 
the current experimental results do not allow conclusive con- 
firmation of the presence and nature of the oscillating struc- 
tures, the isolation of the contribution of /7-H2 in the cross- 
beam experiment or the complementary measurement of dif- 
ferential cross-sections could help to asess their existence and 
origin [73]. If they were resonances, they would be very sen- 
sitive to the potential energy surface and they would give us 
important information about it. Note that due to the high de- 
gree of averaging into orientations, energies and directions of 
approach, direct observation of resonant structures in reaction 
cross-sections at thermal energies h as p roven effusive lf74ll . Up 
to the date it was only successfu HI75ll in the F+HD system, 
and thanks to the concurrence of a set of favorable dynami- 
cal conditions. The current experimental access to the low- 
temperature regime, charaterized by small partial-wave aver- 
aging, may lead to a different scenario. 

Finally, our work reinforces the idea that an accurate treat- 
ment of non-adiabatic dynamics at LR for asymptotically de- 
generate systems seems fundamental at low energies. Previ- 
ous studies at thermal energies stressed the redistribution of 
asymptotic flux onto the different PES which may result from 
such degeneracy while the reagents approach! 64]. The role of 
such non-adiabatic connections at low energies has not been 
explored yet. A better description of the title process includ- 
ing non-adiabatic effects is in progress. 
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VII. APPENDIX : LONG-RANGE ELECTROSTATIC AND 
DISPERSION ENERGIES 

We consider below the case of LR interactions between 
an open-shell atom, A, and a closed-shell diatom, B. To ac- 
count for the n-fold degeneracy of the atom, we define a set 
of asymptotically degenerate diabatic states, taken as a prod- 
uct of the atomic |L/l) (A) and diatomic |0) (B) unperturbed elec- 
tronic states. The quantum number A - -L, . . . , L is the pro- 
jection of the atomic orbital angular momentum L along the 
Body-Fixed (BF) z-axis, chosen along the intermolecular vec- 
tor R, and the projection A — of the diatomic orbital angular 
momentum refers to the diatom axis r, chosen to be in the BF 
xz-plane. 

The multipolar expansion of the first-order, n~?, and 
second-order, H , , perturbation operators depend on multi- 
pole moments Q, and Q ( , B) that are defined in frames cen- 
tered on A and B, with axes parallel to a global Space-Fixed 
(SF) frame. To treat the atom-diatom case, we introduce mul- 
tipole moments Qj^ 2 mat are defined relative to a local refer- 
ence frame on B, 

where Clg collects the polar angles of the diatom axis r relative 
to the SF frame, and D^ 2 „ 2 (Ob) is a Wigner rotation matrix. 
By rotating the SF frame such as the z-axis coincides with R, 
and r lies in the xz-plane, the Wigner rotation matrix reduces 

VI AT* I 

to Z)^ JjUj (Ob) = «m2ju 2 (0)> w i tn @ being the angle between R 
and r. 

Following previous works on open-shell interacting sys- 
tems lE^Hl, the matrix elements of in the \LA) W \Q) <B) 
basis set describe the electrostatic energies and, according to 
the chosen orientation for the global frame, they can be writ- 
ten as: 

V^{R, 9) = J] Vf x i_ m C l2 ,- m (0, 0) (13) 

l\ litn 

where the angular functions Ci yln (6, (p) are normalized spheri- 
cal harmonics, and the electrostatic interaction coefficients are 
given by: 

(2L 
(hml 2 -m\L 12 0) (14) 

X(IA\Qfl\IA') <0|jSg|0> 
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where Lyi = h + h- From the Wigner-Eckart theo- 
rem, it follows that the multipole moments matrix elements 
{LA\Q { ^JLA') and <0|g^J0> vanish unless m = A - A' and 
fi2 = 0, respectively. For the sake of clarity, the electrostatic 

will be substituted hereafter by 



interaction coefficient Vfj 

hhm 



the label Cf^ +l , since it is commonly used in the context of 
LR interactions. 

0)^° basis set de- 



The matrix elements of Hj in the |L/l) (A) 



scribe the dispersion energies, and write as 



ish unless M - A - A' and Q2 = 0, respectively. Further- 
more, in the particular case of dipole-dipole polarizabilities 
(h = 1% — 1) an d a diatomic state {\A\ = \A'\), the matrix 
elements u> a„ , „ n vanish unless Q2 = A - A' and ki is 
even. Again, for the sake of clarity, we substitute hereafter 
the dispersion interaction coefficient Vf^, ir . k M by the more 

commonly used label C$ h+r , +JJci . 



(15) 

where the dispersion interaction coefficients are given by: 
^ 2W = Z4w ^^-Ml^^^ (16) 

k,k 

The angular coupling coefficient ff^,y, is defined as in Eq. 21 
of Ref. [50], and the quantity K^^n ^ is a coupled form 
of Casimir-Polder integral: 

1 r°° 

(17) 

From the Wigner-Eckart theorem, it follows that the cou- 
pled dynamic polarizabilities M ' o^i>) k ^ M and ° a (*/'),t 2 g 2 van " 



In present work, we have considered the quadrupole- 
quadrupole interactions (l\ — I2 = 2) between an open-shell 
atom A and a closed-shell diatom B. Hence, the quadrupole 
moments (LA\Q ( ^_ A ,^\LA'} and (0|<2 2 o'|0) are me necessary 
ingredients to compute the LR electrostatic coefficient Cj 
of Eq. [14] Besides, we have considered the dipole-induced 
dipole-induced interactions (l\ = l' x = li = Zl = 1) b e_ 
tween the two species. In such a case, the LR dispersion co- 
efficient CjH of Eq. [16] vanish unless k\ , &2 and k take even 
values, and the necessary ingredients are the coupled dipole- 
dipole dynamic polarizabilities M ' ahl^a-x) an d ° a (fi)i 2 o' 
with k\ — 0, 2 and £2 = 0, 2. Notice that the matrix elements 
are here defined in terms of spherical components of the Qi m 
multipole moment and ai,„i> m > dynamic polarizability opera- 
tors, together with a basis set of (complex) signed-/! electronic 
states. Nonetheless, we can benefit from the Wigner-Eckart 
theorem to consider only the diagonal matrix elements, which 
have one-to-one correspondence with the quantities derived 
from quantum chemistry calculations. 
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